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ABSTRACT 

The merger of two white dwarfs may be preceded by the ejection of some mass in "tidal tails" , creating 
a circumstellar medium around the system. We consider the variety of observational signatures from 
this material, which depend on the lag time between the start of the merger and the ultimate explosion 
(assuming one occurs) of the system in a Type la supernova. If the time lag is fairly short, the 
interaction of the supernova ejecta with the tails could lead to detectable shock emission at radio, 
optical, and/or x-ray wavelengths. At somewhat later times, the tails produce relatively broad NalD 
absorption lines with velocity widths of order the white dwarf escape speed (~ 1000 km s~ ). That 
none of these signatures have been detected in normal SNe la constrains the lag time to be either very 
short (< 100 s) or fairly long (> 100 yr). If the tails have expanded and cooled over timescales ~ 10 4 yr, 
they could be observable through narrow NalD and Call H&K absorption lines in the spectra, which 
are seen in some fraction of SNe la. Using a combination of 3D and ID hydrodynamical codes, we 
model the mass- loss from tidal interactions in binary systems, and the subsequent interactions with 
the interstellar medium, which produce a slow-moving, dense shell of gas. We synthesize NalD line 
profiles by ray-casting through this shell, and show that in some circumstances tidal tails could be 
responsible for narrow absorptions similar to those observed. 
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1. INTRODUCTION 

Type la Supernovae (SNe la) are important cosmolog- 
ical tools due to their remarka bly standard light curves 
(|Pskovskiil 1197ft iPhillipj 11993ft which lend to their use 
as sta ndard candles (|Colgatdll979l : iBranch fc Tammannl 
1992). They are thought to be the result of a thermonu- 
clear runaway in a degenerate white dwarf star, produc- 
ing large amounts (~ 0.6 M ) of radioactive 56 Ni, how- 
ever the precise mechanism for producing this detonation 
remains elusive. 

Currently, there are two major candidates for SNela 
progenitor systems. In the single degenerate sce- 
nario, the system consists of a white dwarf that ac- 
cretes mass f rom a main- sequence or evolved com- 
panion star (IWhelan fc IberJ 119731 : iNomotol I1982U 
iHillebrandt fc Niemeverl l2000f ). Under certain con- 
straints on the mass-accretion rate, the white dwarf 
mass can near the Chandraskehar limit and a th ermonu- 
clear r unaway will be ignited near the center (jNomotol 
Il982bl ial). In the double degenerate merger scenario, two 
white dwarfs in a close binary inspiral and coalesce. If 
the merger is violent enough, a detonation could be trig- 
gered on the dynamical timescale (~ 10 — 100 sees) , 
promptly exploding the sy stem (jPakmor et al.l I2010L 
120111 : iGuillochon et al.ll2010D . Absent a detonation, the 
resulting configuration is a more massive, degenerate 
remn ant surrounded by an accretion disk (jYoon et al.l 
120071 ). The disk will evolve viscously on timescales 10 4 — 
10 8 s, heating the shear l ayers and driving the remnant 
to a more spherical state (Shen et al. 2012; Sch wab et al.l 



120121 : Uiet aTll2013[ ). The envelope will then evolve ther- 
mally on timescales approaching 10 3 — 10 4 years, while 
the core thermal timescale is ~ 10 5 years. If the sec- 
ondary white dwarf is composed of helium, stable shell 
burning may ensue on similarly long timescales. 

The final outcome of such a merged remnant is stil l 
debated. Pioneering calculations (Saio fc Nomotoll2004f ) 
indicated that system likely evolved to an off-center, non- 
violent ignition, with the end result being collapse to a 
neutron star. Subsequent studies have challenged this 
conclusion, suggesting that the remnant may instead 
evolve to central i gnition (and a SN la) on a viscous or 
thermal timescale (Yoon et al.l 120071 : Ivan Kerkwijk et all 
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120101 : Ui et al.ll2013f ). Other detailed studies of the post- 
merger evolution, however, suggest that off-center ig- 
nition and colla pse is indeed the more likely re sult in 
most cases (jShen et all 120121 : iSchwab et al.ll2012h . The 
remnant evolution apparently depends on the binary 
mass ratio and the initial synchronization conditions 
(|Zhu et al.ll2012l ). and may be sensitive to the treatment 
of e.g., mass loss - hence, thermonuclear explosions may 
be allowed in certain limited regions of parameter space. 
In addition, entirely different ways of triggering a deto- 
nation in a double d egenerate system may be relevant. 
iKatz fc Pond (|2012[ ) suggest that if the binary is in a 
triple system, Kozai resonance effects can lead to multi- 
ple close passages (with period of ~ 10 4 y rs) culminating 
in a near head on collision of the stars (iRosswog et all 
120091 : iRaskin et al"1l2QQ9L 120101 : lHawlev et al.ll2012l ). 

Given the difficulty of the theory, it would be valu- 
able to identify empirical signatures of double degenerate 
mergers. One of the more dramatic consequences of the 
dynamics of compact object mergers is the tidal strip- 
ping and ejection of mass which occurs just prior to coa- 
lescence. These "tidal tail" ejections are a robust feature 
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of merger c alculations and c onsist of 10 ~ 4 — 10 ~ 2 M 
of material (jDan et al.l l2011f ) moving at the escape ve- 
locity, ~ 10 8 cm s _1 , and concentrated in the equatorial 
regions. Though the total mass expelled is a small frac- 
tion of the system mass, it is still enough to produce a 
relatively dense surrounding medium. If the circumstel- 
lar material (CSM) is relatively nearby (r < 10 16 cm) 
when a SN la explodes, the interaction of the supernova 
ejecta could lead to detectable shock emission at radio, 
optical, and/or x-ray wavelengths. If the CSM has in- 
stead expanded to larger radii and cooled, it may still be 
observable through narrow absorption lines in the spec- 
tra, in particular the resonance doublet lines of NalD and 
Call H&K. 

Recently, spectral observations of a handful of 
SNela have discovered such narrow CSM absorp- 
tions a t velocities (~ 1 — 100 km s~ ) (|Patat et al.l 
20071: iB orkows ki et all 120091: ISimon et aD 120091 : 
Sternberg etaJ] 120111 : iFolev et al.l 12012( 1. In a few 



cases, the strength of the lines are observed to vary 
with time, suggesting that the CSM is relatively near 
to the supernova. In most other cases, the lines do not 
vary with time, but are observed to be prefe rentially 
blue-shifted for SNe la ([Sternberg et al.l 1201 lh . which 
suggests that they are statistically related to mass-loss 
from the SNela progenitor system. Typically, single 
degenerate systems are assumed to be the culprit for this 
CSM, the result of wind mass loss from a red giant or 
asymptotic giant branch companion, perhaps shaped by 
novae eruptions f rom the whi t e dwa rf ([Moore & Bildstenl 
12012( 1. Recently, IShen et al.l ([20131 ) have suggested that 
for certain classes of double degenerates (C/O + He 
WDs) nova eruptions in accreted surface layers could 
produce the necessary absorber material. We show here 
that the tidal tails from double degenerate mergers, 
when combined with interactions with the interstellar 
medium (ISM), can also result in a dense cloud of 
blue-shifted, absorber material moving at a rage of 
velocities, assuming that the time-scale between the 
mass ejection and the SN la explosion is long. 

In this paper, we explore the signatures of CSM from 
tidal tails in double degenerate mergers. In §2, we discuss 
the combination of 3D and ID numerical approaches we 
use for modeling the white dwarf merger as well as the 
subsequent interaction of ejected material with the ISM. 
We also compare the results of these numerical studies to 
analytical estimates of this interaction in §2.1. In §3, we 
estimate the optical depth of the absorbing material and 
use a spectral synthesis code to produce an absorption 
profile of the NalD line. Finally, in §4, we summarize our 
conclusions and discuss potential avenues for followup. 

2. TIDAL TAIL EJECTA 

To determine the mass-loss during a double degener- 
ate white dwarf merger, we use a 3D smoot hed particle 
hydro dynamics (SPH) code called SNSPH ([Frver et al.l 
|2006[ ). SPH codes are lagrangian by design, and so nat- 
urally conserve angular momentum. This feature is cru- 
cial for estimating the fraction of the donor star that at- 
tains escape velocit y. Following the procedure laid out in 
iRaskin et al.l (|2012[ ) for generating accurate initial condi- 
tions for a white dwarf merger, we simulate a merger pair 
consisting of a 1.20 M primary and a 0.64 M secondary 
with 5 x 10 5 equal mass particles throughout each star. 



This mass pair is disparate enough for the secondary to 
experience strong tidal forces, while still close enough in 
mass for the merger to exhibit catastrophic disruption of 
the smaller mass star. 

As the secondary star disrupts, material flows through 
the L2 point and escapes the system. Figure [1] illus- 
trates the geometry of the simulation shortly after the 
complete breakup of the 0.64 M donor star into a disk. 
Roughly 2.4 x 10 -3 M of material achieves escape veloc- 
ity (« 2000 km s _ ), forming the end of a long tidal tail. 
The ejecta has an opening angle of ~ 93° in the plane of 
the disk and « 41° in the perpendicular direction, giving 
a covering factor > 25% seen edge on, and ~ 9% for all 
space. 




Figure 1. A density slice snapshot of the 3D SNSPH simulation of 
a white dwarf merger of masses 0.64 Mq + 1.20 Mq . Material from 
the donor star that is no longer gravitationally bound is enhanced 
in brightness and colored blue. The ejecta material is represented 
here by 806 particles in this simulation. 

Simulations of other mass pairs yield very similar tidal 
tail ejections. Table [T]lists the ejecta masses from each of 
our SPH simulations, and as is evident from these results, 
tidal tail ejections massing M e j ~ 10 -3 M are indeed a 
robust feature of double degenerate white dwarf mergers, 
irrespective of the white dwarf masses. Our analysis of 
the subsequent evolution of the tidal tails will be confined 
to our fiducial case of 0.64 M + 1.20 M . 



Table 1 

Ejected masses, M e j, from various SPH simulations of white 
dwarf mergers of masses Mi and M<i . 



Mi/M 


M 2 /M 


M ei /M 


0.64 


0.96 


2.0 x 10- 3 


0.64 


1.06 


2.3 x 10~ 3 


0.64 


1.20 


2.4 x 10~ 3 


0.84 


0.96 


1.4 x 10- 3 


1.06 


1.06 


4.7 x 10- 3 


1.06 


1.20 


3.3 x 10- 3 



Shortly after its ejection, the tidal tail expands bal- 
listically within the nearly central potential field. It is 
expected that this material (at roughly 0.6 g cm -3 in 
the snapshot represented in Figure [1]) should drift rela- 
tively unimpeded for ~ 100 years before the swept up 
ISM mass becomes non- negligible. During that time, it 
will expand considerably and drop in density by many 
orders of magnitude. Since an SPH simulation of this 
expansion would be costly, we instead evolve this ma- 
terial in a simple n-body code with self-gravity until it 
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reaches homology. The merger remnant (all bound par- 
ticles) is treated as a point mass for this phase of the 
calculation. As the SPH particles were of roughly equal 
mass (ra), the density at any spatial point r in the n-body 
calculation can be approximated by m f -3 , where f is the 
average interparticle spacing around r. Once the mate- 
rial is expanding homologously, rescaling the domain is 
trivial. 

When the material has drifted far enough to have 
reached a region where ISM interaction becomes relevant 
(~5x 10 17 cm), we map the resultant distribution into 
a ID spherical lagrangian hydrodynamics code with a 
gamma-law equation of state. Figure [2] shows the initial 
conditions for this phase of the calculation. Although the 
tails are clearly aspherical, the interaction and sweeping 
up of the symmetric ISM likely causes the distribution to 
spread and approach a more spherical configuration, in 
which case a ID description is likely not too inaccurate. 
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Figure 2. The initial conditions of the ID lagrangian simulation 
for M e j = 2.7 x 10 -3 Mq . Densities here are plotted as proton 
masses per cm 3 (1.67 x 10 -24 g cm -3 ) to better illustrate the 
density contrast of the ejecta with the unshocked ISM. The ejected 
material has expanded homologously to a density of ~ 3.7 x 10 -22 
g cm -3 and is just beginning to contact the ISM at v r ~ 2000 
km s — 1 . 

The typical density of the ISM is uncertain, however 
estimates of the ambie nt medium dens i ty aro und Tycho's 
supernova remnant by iBadenes et al.l (|2006f ) suggest an 
ISM density of between 0.2 x 10" 24 g cm -3 and 5 x 10" 24 
g cm -3 , with the most probable density at 2 x 10 -24 g 
cm -3 . This is consistent with roughly 1 proton mass per 
cm 3 , and for our ID, late-phase simulation, we use this 
conservative estimate, /?i sm = 1.67 x 10 -24 g cm -3 . 

Since radiative cooling of the shocked ejecta material 
is important over the long time-scales considered here, 
we employ an analytical cooling function of the form 



A 1.35 x 10- 16 , A in 27 r- 

= = + 1.4 x 10~ 27 g^VT 

n e rii T 



(1) 



for temperatures above 10 5,5 K. In this expression, 
n e and rii are the number densities of free electrons 
and of the ions, respectively, and g^ is the Kramers- 
Gaunt factor which can safely be set to 1.0. The 
second te rm in this expression is based on free-free 
emission (jRvbicki fc Lightmanl Il986[ ) for temperatures 
above ~ 10 75 , and the first term approximates the 
theoretical cooling curves from emission lines found in 



iSutherland fc Dopital (|1993f ) and IWiersma et al.l ([20091 ). 
For simplicity we take A/(n e rii) = 1.35 x 10 -21,5 
erg cm 3 s _1 for 10 4 < T < 10 5 ' 5 K. This approximates 
the collision ally ionized equ il ibrium cooling rates in this 
range from IWiersma et al.l (J2009) fairly well, however 
sodium is expected to remain ionized at these temper- 
atures. Since cooling rates below 10 4 K are somewhat 
undetermined, and since the exact functional form of the 
cooling function at this low temperature does not impact 
the structure of the shell, we use the same constant rate 
down to 10 3 K. 

Our simplified cooling function captures most of the 
relevant cooling physics for the shell formation with a 
slight decline in cooling efficiency from 10 5,5 K out to 
« 10 T K and a steady rise in efficiency thereafter, while 
also being efficient enough at low temperatures to al- 
low for neutral sodium formation. The extent to which 
we overestimate the neutral sodium fraction is only im- 
portant insofar as our inferred optical depths will be too 
great, and determining precisely what the sodium ioniza- 
tion fraction should be is out of the scope of this work. 

When the ejecta first encounters the ISM, it is shock- 
heated and begins to slow as it sweeps up ISM gas. Be- 
hind the shock, the mixed ISM and ejecta material ini- 
tially cannot cool, and the hy drodynamic evolution is 
essentially a Sedov blast wave ()Sedovl[l959[ ). However, 
as the post-shock gas expands, the temperature drops 
sufficiently for collisional ionization to begin to cool the 
gas, and after w 3 x 10 3 years, a thin, dense shell forms 
near the shock front. This shell has a mass of ~ 1.5 M 
and a velocity < 40 km s _1 . Figure [3] shows the growth 
of this dense shell. 
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Figure 3. A later stage of the ID lagrangian simulation for M e j = 
2.7 x 1CT 3 M . The shocked ISM gas has cooled to T < 10 3 K 
and formed a dense shell moving at v r < 40 km s — 1 . Densities 
here are plotted as proton masses per cm 3 (1.67 x 10 -24 g cm -3 ) 
to better illustrate the density enhancement of the shell. 



2.1. Analytical Comparison 

To compare these results with analytical expectations, 
we consider the three phases for the evolution of the ex- 
panding gas. The first is a simple ballistic phase where 
r e j = v e jt. The duration of this phase is roughly the time 
before the ejecta sweeps up its own mass of ISM gas, 
and lasts approximately 200 years. At this point, the 
gas will enter an energy conserving (Sedov) phase where 
the temperatures of the shocked shell are high and ra- 
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diative losses are negligible. Evenutally, when the gas 
has cooled sufficiently to become radiative, it will enter 
a momentum conserving phase. 

During the energy conserving phase of the interaction, 
when the adiabatic evolution is t hat of a Sedov blast 
wave (jSedovl 119591 : lChevalierl[l982h . the location of the 
contact discontinuity is, for a uniform ISM, roughly of 
the form 



E e] t 



1/5 



0.87 



P'lsm 



1/5 



(2) 



where E e - } is the kinetic energy of the ejecta material. 
The shock velocity then evolves as 



^Mfvfpr^r^. 



(3) 



The jump conditions for a strong shock give a post-shock 
pressure of the form 



Ipi* 



2 7+1 



^^ism^-^ ,s« 



(4) 



Here, 7 is the usual adiabatic index, and the post-shock 
temperature for an ideal gas is then 



T s « 1.2 x 10 7 /i^ 8 K 



(5) 



where v s ^ = v s /10 8 cm s 1 . The timescale for this post- 
shock gas to cool is 



t'COOl 



nkT s 

(7-1) A' 



(6) 



Employing an approximation for the cooling function 
which applies in the range 10 5 K< T < 10 7,3 K from 
lDrainel(|20TTl ). A /(n.nv) w 1.1 xlO" 22 ^ -0 - 7 erg cm 3 s"\ 
and combining the previous relations gives for a nominal 
ISM number density of 1 cm -3 



tcool^^OT^S 



^2.3xl0 34 M^3 68 ^ 3 8 6 t- 2 - 04 s, 



(?) 



where M_ 3 = M ej /1O" 3 M and v e j,8 = ^ej/1000 
km s _1 . The time to cool (t coo i ~ i) is then roughly 
10 4 years. 

At this stage, the gas becomes radiative and energy is 
no longer conserved. The shell will reach a velocity Vf 
given by conservation of momentum 



M ej ^ ej 






shVf, 



(8) 



where M s h is the final mass of the shell, which for the 
case Vf <C v e j will be dominated by swept up ISM gas. 
Assuming the ISM is of constant number density, ni sm , 
the radius at which the shell slows to a value Vf is then 



?"sh = 



M e 



47r /ira p n ism Vf 



1/3 



(9) 



!M l/3 1/3 -1/3 -1/3 
/lvl -3 V ej,8 V f,6 r 



1/3 
ism P C 5 



where Vf^ = Vf/10 km s and n* sm is the dimensionless 
quantity ni sm /cm -3 . For our purposes here, these quan- 
tities are all of order unity, and the predicted location 



of the shell is ~ 1.1 pc. Taking the time derivative of 
equation (|9]) gives the drift time 



1 ( 3 MejVej 
Mrift — T 



4 \ 47r fjLm p ni S 



1/3 



-4/3 



v f 



(10) 



where V{ is the velocity at the beginning of the radiative 
phase. This drift time is also ~ 10 4 years, and so the 
total time from the ejection of the tidal tail to the inter- 
action with the ISM and subsequent slowing to a cooled, 
radiative shell moving at ~ 40 km s _1 is w 2 x 10 4 years. 
This compares well with the simulated result of 1.9 x 10 4 
years. 

Figure [H diagrams the evolution of the shock through 
the three phases outlined here, comparing the analytical 
estimates to the simulation results. The analytical esti- 
mates are only relevant after an initial drift period, dur- 
ing which the ejecta mass is much greater than the swept 
up ISM mass. The shock radius during the energy con- 
serving phase evolves as a 2/5ths power-law, and during 
the cooling phase, as a l/4th power law. At this stage, 
the shell velocity in the simulation becomes difficult to 
measure accurately as without clumping and self-gravity, 
it begins to expand adiabatically into the surrounding 
medium. 
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Figure 4. The evolution of the ejecta mass interacting with the 
ISM, forming an absorbing shell. Both the results of a ID la- 
grangian simulation and the analytical estimates at At e x > 100 
yr are shown here. The phases illustrated here are separate from 
those discussed in 52 and are discussed in more detail in 83. 



3. OBSERVATIONAL SIGNATURES 

The potentially detectable signatures of the tidal tails 
depend primarily on the lag time, At ex , between when 
the tails are ejected and when the merged system ex- 
plodes as a SNe la (assuming that it does explode). The 
value of At ex is unclear given our incomplete understand- 
ing of the post-merger evolution, but obviously interest- 
ing physical timescales include the dynamical (At ex ~ 
10 2 - 10 3 sees), the viscous (At ex ~ 10 4 — 10 8 sees), and 
the thermal (^ 1Q 4 — 1Q 5 years) time s. For resonantly in- 
duced collisions (jKatz fc Dongll2Q12[ ) a relevant timescale 
is the period of close approaches (~ 10 4 years). Without 
adopting any specific theoretical paradigm, we consider 
the possible observable signatures as a function of At ex . 

3.1. At ex < 100 s: shock breakout and cooling 
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If the system detonates promptly, on the dynami- 
cal timescale, the tidal tails will have only had time 
to reach radii of r < 10 10 cm before they are im- 
pacted by the supernova ejecta. At these times, the 
CSM is optically thick, and the interaction will drive 
a radiative shock through it. Emergence of this shock 
should produce a brief x-ray pul se similar to shock break- 
out in core collap se su pernovae (ICh evalier & Klein 19791 : 
llmshennik et all 119811 : iHoflich fc^chae fer 20091) . Based 
on scaling relations derived fo r x-ray bursts in SNe la 
(jPiro et al.l 120101 ; iKasenl l2010[ ) we anticipate bolomet- 
ric luminosities of order 10 44 ergs s _1 , peaking in the 
1 — 10 keV range and lasting roughly a light crossing 
time, t ~ 10 _2 At ex ~ 1 sec. 

Subsequent cooling of the shock heated layers will con- 
tinue to produce UV/optical emission for some time after 
breakout, although this will be relatively dim and diffi- 
cult to detect unless the supernovae is caught extremely 
early. The cooling emission will also be anisotropic, and 
may only be prominent for viewing angles of the half of 
the system where the ti dal tail is ejected. The early ob- 
serva tions of SN 2011fe (jNugent et al.ll20 lit [Bloom et al.l 
2012) limit the emitting radius to be < 1.5 x 10 9 cm, 
which corresponds to a lag time At ex ~ 10 s. This con- 
straint assumes that the explosion time has been accu- 
rately determined b y extrapolating the e arly light curve, 
which is not trivial (jPiro fc N akar 2012). More detailed 
modeling is required, but it is likely that the early obser- 
vations of SN 201 lfe and other SNe la (e.g., SN 2009gi, 
Foley et al. 2012) are sensitive to tidal tail interaction 
and hence constrain WD detonations from violent merg- 
ers. 

3.2. 100 s < At ex ;$ 10 6 s: early UV/optical emission 

For this range of lag times (comparable to the vis- 
cous timescale of the remnant disk) the CSM is confined 
to radii r < 10 13 cm and will be impacted by the su- 
pernova ejecta within a few hours after the explosion. 
The resulting shock breakout burst will be longer lasting 
(^minutes) and fairly soft (~ 100 eV). Following this, 
optical/UV cooling emission from the expanding, shock- 
heated layers should be quite bright, Lbol ~ 10 42 — 10 43 
erg s _1 . Statistical studies of early SN la light curves 
have strongly ruled out shock emission o f this magni- 
tude i n the first few days after explosions (Bi anco et al.l 
l2Qllfc lHavden et"a l. 2010). However, because the total 
mass of the tidal tails is relatively small, the duration of 
that emission should be fairly short - of order the dif- 
fusio n time through the shocked region, or about half a 
day (jPiro fc Nakarll2012[ ). Emission this early may have 
been misse d in the statistical s tudies. For the ca se of 
SN 2011fe (jNugent et aLll20Tfl : iBloom et alJl2012h . the 
early light curve data clearly rules out tidal tail interac- 
tion at these radii, unless one believes that the explosion 
time had been underestimated by ~ 0.5 days. 

3.3. 10 6 s < At ex < 10 8 s: interaction 

For these lag times, the CSM is spread over radii 
r < io 14 — 10 16 cm, and the supernova ejecta may ex- 
perience interaction over timescales of weeks. One ex- 
pects a range of signatures, in particular the shocked gas 
should produce radio and x-ray radiation due to syn- 
chrotron and inverse compton scattering, while the pho- 



toionization of the unshocked CSM may produce optical 
emission lines. Relatively speaking, the density of the 
tidal tail CSM at these phases is high, with an effective 
mass loss rate in the range M w M e j(v e j/v w )/At ex w 
10 -2 — 10 -5 M yr _1 depending on the tidal tail mass 
and the wind velocity, v w . This is significantly greater 
than the CSM density of most companion star winds in 
the single degenerate scenario. 

A variety of SN la observations at x-ray/optical/radio 
wavelengths have constrained the mass loss rate to be 
5- 10~ 6 — 10~ 5 M^ yr -1 for several individual events 
(iHughes et al.ll2007l: i Mattila et al.ll2005l: iCumming et al.l 
119961 : iPanagia et al.ll2006l : iRussell fc Immlerl 120 12( 1. The 
observations of SN 201 lfe are the most constraining, 
limiting M < 10~ 8 M^ yr" 1 (iHoresh et a p 120121 : 
iMargutti et~aLll2Q12l : iChomiuk etHI I2012D . These re- 
sults apparently exclude white dwarf merger models with 
lag times in this range. 



3.4. 10 8 s < At e 



< 100 yrs: "broad" absorption 
features 



Over these lag times, the bulk of the CSM has ex- 
panded ballistically to r ~ 10 17 cm and it will take years 
for the supernova ejecta to interact with it. However, one 
may still be able to see the tidal tails in absorption. The 
tails will have cooled adiabatically and are likely neutral 
with Nal columns of order 10 13 - 10 15 cm -2 . The NalD 
lines will then be optically thick and will produce ab- 
sorptions for those viewing angles where we are looking 
through the tails, or about 10% of the time. The width 
of these CSM absorptions will be of order the escape 
velocity, ~ 10 8 cm s _1 , which is narrower than typical 
supernova lines, but much broader than typical CSM ab- 
sorptions. There is therefore some hope of uniquely iden- 
tifying lines associated with tidal tail ejections. NalD ab- 
sorption features of this width have not been identified 
in any SNe la, despite the enormous number of optical 
spectra obtained. This would apparently rule out white 
dwarf merger remnants that explode with these lag time 
scales. 

3.5. 100 yrs < At ex ^ 10 3 yrs: no signature 

After roughly 100 years, the expanding tails have swept 
up enough ISM material that they have been entirely 
shocked and are in the Sedov phase. The hot shell of gas 
will be ionized, and optically thin to all relevant lines. 
One therefore does not expect any prominent observa- 
tional signatures again until the shell has begun to ra- 
diatively cool. 

3.6. 10 3 yrs < At ex < 10 5 yrs: narrow absorption lines 

After 10 3 yrs, the tidal tail ejecta have slowed to 
the point that radiative cooling becomes important and 
the CSM forms a dense shell moving at velocities < 
100 km s~ . The interaction with the CSM likely causes 
the tails to spread and take on a more spherical shell 
configuration. Absorption lines may then be seen from 
the cool shell, at least in a statistical way. We can cal- 
culate the optical depth of the NalD line (at line center) 
through a shell of density n s h and thickness Ar as 



r = n sh x Na ia Ar, 



(11) 
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where XNai is the fraction of sodium in the neutral state 
and (Jo is the cross-section, per atom, of the NalD line at 
line center. Assuming that the width of the line is set by 
Doppler broadening at a velocity Av ~ 10 km s~ , this 
cross-section is 



Tre 2 A 

&0 = T— Josc^Nal^Na 

m e c Av 
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Av 6 1 x Na i cm 2 , 



(12) 

where we take the number abundance of sodium to be 
solar, ^Na = 2 x 10 -6 . The neutral fraction of sodium, 
XNai is difficult to estimate, since it likely requires cooling 
below the recombination temperature of hydrogen. In 
the absence of detailed calculations, we assume here that 
the neutral fraction is of order unity. 

For a thin shell of thickness Ar, the shell density is 
p s h = M s h/(47rr 2 h Ar) and the line optical depth is in- 
dependent of the shell thickness. Using equations ([TT]L 
(fl2|h and (|9j), and given that the swept up ISM mass 
dominates the mass of the shell, we can then write the 
optical depth of the Nal lines in terms of the density of 
the ISM and the final velocity of the shell as 



{ lQ7r 2 M e] nl m (v e] /v f ) \ 
9[im p J 



1/3 



X lO-^A^^Nal 



:1.39 X Na I 



v /,6 



1/3 



(13) 



For the values used in this paper (n{ sm = 1 and Vf^ = 
4), this reduces to r ~ 0.87. That the optical depth 
depends only on the ISM density and the final velocity of 
the shell should not be too surprising as an undisturbed 
ISM will reach an optical depth of ~ 1 for sodium at a 
radius where Vf$ ~ 5. Without the moving shell, the 
Na absorption would appear as narrow ISM lines (not 
blue-shifted), as is seen for other transient types. Put 
another way, if the source of the blue-shifted absorption 
were from a shell of shock-cooled ISM gas (as is the case 
here since M e j/M s h << 1), requiring an optical depth 
~ 1 sets a radius (similar to the same requirement for 
unshocked gas, modulo a factor of 3) at which Vf$ « 5 
naturally, provided M e] is sufficiently large for equation 
(0) to hold. 

In order to synthesize an absorption line profile for 
our simulations, we simply integrate equation ([TT]) nu- 
merically for every wavelength A and use a gaussian 
dispersion about the line center Ao for thermal broad- 
ening, A A = (Xo/c)(^/kT//jLmp). We take XNai, the 
neutral fraction of sodium, to be 1.0 everywhere where 
T < 5 x 10 3 K. Figure [5] shows the evolution of the neu- 
tral sodium absorption line through various phases of our 
simulation. 

In a simple model, such as this, the shell cannot persist 
for long as it is at a higher pressure than its surround- 
ings. Without clumping due to self-gravity, the shell sim- 
ply begins to expand into the surrounding medium and 
the subsequent drop in density precludes the synthesis 
of any strong absorption lines at much later times. It 
is expected that clumping from Raleigh- Taylor instabili- 
ties not captured in a ID simulation should also enhance 
absorption at later times. These factors work to con- 
strain the available time in our simulation for sodium 
absorption. However, we expect a multi-D approach in- 
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Figure 5. Sodium absorption profiles for a range of times from 
the early phase of the formation of the shell, when the shell has 
not had sufficient time to cool in order for neutral sodium to form, 
to the narrow absorption line phase, to the late-time disruption of 
the shell due to adiabatic expansion. At very late times, the zero 
velocity absorption material has all been swept up in the finite- 
sized simulation box. 

corporating self-gravity should, for the most part, vali- 
date our results and extend the available absorption time. 
If we compare our synthesized line profiles from times be- 
fore adiabati c expansion to observed absorption lines for 
SN 2008fp in lSternberg et al.l ([201 lh . we find a good cor- 
respondence for the absorption profile of the blue-shifted 
sodium lines. 

4. DISCUSSION 

We have explored the observational consequences of 
the tidal tails believed to be ejected in SNe la resulting 
by white dwarf mergers. The tails produce several ob- 
servational signatures that may constrain the lag time 
between the merger and the subsequent explosion of the 
remnant (assuming an explosion does occur). 

If the lag time is very short (A£ ex < 100 s), the inter- 
action of the supernova ejecta with the tidal tail ejecta 
should produce brief x-ray pulses like those observed dur- 
ing the shock breakout phase of core collapse supernovae. 
If the lag time is slightly longer (At ex ~ 10 4 s), the inter- 
action shock results in early (~day) UV and optical emis- 
sion. For lag times > 10 8 s, the tidal tail ejecta material 
is too distant to be impacted by the supernova ejecta in a 
meaningfully short time, but the tail material may still be 
seen in absorption features that are broader than typical 
CSM absorption lines. The lack of observed shock emis- 
sion or broader CSM absorption features in the SN la 
sample would appear to exclude mergers exploding with 
lag times in the range 10 4 s < At ex < 100 yr, which 
places important constraints on models presumed to ex- 
plode on a viscous timescale. 

By combining the results of a 3D hydrodynamics simu- 
lation with simple ID models, we have shown that given 
time to interact with the ISM (A£ ex > 100 yr), these mass 
ejections produce Sedov shock waves in the surrounding 
medium which then evolve into thin, over-dense shells af- 
ter radiatively cooling (A£ ex > 10 4 yr). The covering fac- 
tor for the ejecta is < 10%, however, interaction with the 
ISM will likely spread the tail material and increase the 
covering fraction to something closer to ~ 50%. Multi- 
dimensional models will be needed to determine the ge- 
ometry and the sub-structure that develops due to hy- 
drodynamical instabilities. 



Tidal Tail Mass Ejections 



Depending on the exact value of At ex , the resulting 
shock front can have a range of velocities from 10-100 
km s~ , reproducing the observed range of blue-shifted 
sodium absorption lines in SNela. The timescale for the 
evolution of this absorption material into a shell is also 
comparable to the timescale for the thermal evolution or 
helium shell burning time of white dwarf merger rem- 
nants, ~ 10 4 — 10 5 years. The optical depths of these 
shells are of order ~ 1 , independent of many of the sim- 
ulation parameters. Synthesizing absorption lines from 
our hydrodynamical simulati ons yields absorption pro- 
files similar to those found in St ernberg et al.l (|2011f ). 

Although we have focused here on CSM due to tidal 
tail ejection, there may be other means of creating 
mass outflow in merging double-degenerate systems. 
These include mass outflows during th e rapid accretion 
that immediately prec edes the merger ([Dan et al.l 120111 : 
iGuillochon et al.ll2010D, disk w inds arising during the vis- 
cous evolution (|Ji et al.l 12013). or mass ejections in the 
post- main-sequence evo lution of the white dwarf progen- 
itors (Shen et al.ll2013l ). Another interesting case is the 
possible tidal stripping that may occur during close en- 
counters in triple sy stems, d riven by Koza i reson ances 
(jKatz fc Dondl20l2h . While iKatzfc Pond (J20121) only 
consider "clean case" scenarios, there is potentially a 
large parameter space available for Kozai resonances to 
induce tidal stripping at the last closest encounter before 
the two of the stars collide. Further simulation will be 
needed to quantify the mass loss by some of these mech- 
anisms. Tidal tail ejection, however, appears to be a ro- 
bust prediction of white dwarf mergers; searching for the 
resulting signatures should therefore provide interesting 
constraints on the progenitors of SNe la. 
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